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Abstract Lame ordinary differential equation in the Weierstrass's form and Heun 
equation are of Fuchsian types with the three regular and one irregular singularities. 
Lame equation in the Weierstrass's form is derived from Heun equation by changing 
all coefficients y = 8 = £ = A, a = p ~ 2 , j3 = — (a + l),q = —hp~ 2 and independent 
variable.* = sn 2 (z,p). (see (1) in Ref. lfTBI ) 

I consider the power series expansion of Lame function in the Weierstrass's form 
and its integral forms applying three term recurrence formula! D- 1 investigate asymp- 
totic expansions of Lame function for the cases of infinite series and polynomials. I 
will show how the power series expansion of Lame functions in the Weierstrass's 
form can be converted to closed-form integrals for all cases of infinite series and 
polynomial. One interesting observation resulting from the calculations is the fact 
that a 2F1 function recurs in each of sub-integral forms: the first sub-integral form 
contains zero term of A' n s, the second one contains one term of A„'s, the third one 
contains two terms of A„'s, etc. 

This paper is 7th out of 10 in series "Special functions and three term recurrence 
formula (3TRF)". See section 7 for all the papers in the series. Previous paper in series 
deals with the power series expansion and the integral formalism of Lame equation 
in the algebraic form and its asymptotic behavior [19]. The next paper in the series 
describes the generating functions of Lame equation in the Weierstrass's form[21 1. 
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1 Introduction 

In 1837, Gabriel Lame introduced second ordinary differential equation which has 
four regular singular points in the method of separation of variables applied to the 
Laplace equation in elliptic coordinates] 10|. Various authors has called this equation 
as 'Lame equation' or 'ellipsoidal harmonic equation' |[TD . 

Due to its mathematical complexity there is no analytic solution in closed forms 
of Lame function! 11, 12, 13 1. Because its solution, in the algebraic form or in the 
Weierstrass's form, was a form of a power series that is expressed as three term re- 
currence relation lfT2][T3l . In contrast, most of well-known special functions consist of 
two term recursion relation (Hypergeometric, Bessel, Legendre, Kummer functions, 
etc). 

In my previous paper[2|, applying three term recmrence formula. Q), I showed 
the power series expansion in closed forms of Lame function in the algebraic form 
(infinite series and polynomial) including all higher terms of A„'s by applying three 
term recurrence formula. [1"|. I obtained representations in form of contour integrals 
of Lame function in the algebraic form and its asymptotic behavior of it and the 
boundary condition for x. 

In this paper I will show the analytic solution of Lame equation in the Weier- 
strass's form. Its functions in the Weierstrass's form appear as we apply the method 
of separation of variables to Laplaces equation in an ellipsoidal coordinate system 
(Gabriel Lame 1837Q3). 

The Lame equation in Weierstrass's form in which is 

^ = {a(a + l)p 2 sn 2 (z,p)-h}y(z) (1) 

where p, a and h are real parameters such that < p < 1 and a > —j. If we take 
sn 2 (z,p) = % as independent variable, Lame equation becomes 

d 2 y 1/1 1 I \ dy -a{a + l)E,+hp- 2 



d$ 2 ' 2 V4 1 4-i 1 i-p~ 2 )di 1 4£(§-i)(5-p- 2 )' y ^ (2) 

This is an equation of Fuchsian type with the four regular singularities: E, = 0, l,p _2 ,<» 
The first three, namely 0, l,p -2 , have the property that the corresponding exponents 
are 0, A which is the same as the case of Lame equation in the algebraic form. In 
Ref.[2], Lame equation of the algebraic form is 

d2 y , 1 ( 1 , 1 , M d y , -a(a + l)x + q =Q 



dx 2 2\x — a x — b x — c/dx A(x — a)(x — b)(x 

If we compare (f2]i with (0, all coefficients on the above are exactly correspondent to 
the following way. 
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We can another expression of Lame function in the Weierstrass's form by using <(4j 
in Ref.Q. 



2 Power series 

2.1 Polynomial in which makes B n term terminated 

The power series expansion of Lame function in algebraic form for the polynomial 
in which makes B„ term terminated in Ref.f2] is 



y(z) 




"o (-00)^(00+1+ A) 



fO~in 



(1+2 )'o(4 + 2 )'o 

(2 a -fc- c )(/ + |) 2 - a (a + |)(ao + I + |) + ^ (-^^(ao + l + A),, 
(/o + Kl)0o + i + |) (l + |)«o(! + |)/o 

2 ' 2^'oU T jjj), 




x 



where 



(-(*!);„(«! + § +A) io (| + f ),-,(§ + |) ;] 

00 (2 a -fc~ C )(/ + |) 2 - a (ao + |)(ao + | + f) + ^(-ao) iQ (aQ + I + A), 
(''o + 3 + t)('o + | + |) (1 + |)<o(! + |)/o 

ga-fc - § + 1) 2 - a (a, + f + |)(a, + | + i + |) + ^ 

(4 + | + 5 + |)(4 + § + 4 + |) 

{- a k)i k {ttk + k+\ + A),- t (1 + 1 + §),,_, (| + | + 1 
{-a k ) ik _ { {a k +k+\ + X)i k _ i (1 + | + + § + J 

f (-^kK+^+i+^kd+f + lk-.d + f + !)-■„-, ,„] ■] 

^(-^(an+n+I+A^a + l + lkd + f + IX J J 



Z = x — a 

= (a-b)(a-c) ^) 
_ ~(x-a) 
r 1 ~ (a-b)(a-c) 



and 



a = 2(2a ; + i + A) or - 2(2a ; + i + A) - 1 where /, a,- = 0, 1,2, 
a, < a ; - only if i < j where i,j = 0, 1 ,2, ■ ■ ■ 
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Put © in ©-(H). And put c = 1 as X=0 and c = (-p 2 ) 1//4 as A = \ into them. Then, 
we obtain two independent solutions of Lame equation in the Weierstrass's form. 

;Ai = -p 2 sn 2 (z,p);ri = -p 2 sn 4 (z,p) 





(a 








M 




-,>,-!( 








\2~ 


v2~~ 



y=o,i,2,. 



<*> (-«o) !0 (ao + !),- . q 



E 

'0=0 



(4)10(1)10 

-(1+p 



^ (*> + *)(*> + *) 



2^1 (-ao)i (ao + ^) 



'0=0 



4/ 

£ r g -( 1 +P" 2 )'o+2^ (-«b)i (ab + i)fo 

«=2\i =0 (*0 + |)(*0 + |) 



(4)10(1)10 



A (-«i).b(«i+!).b(!)ii(!)ii 



(itk(ik 



*n( e _(1+r2)(,i ^ )2+ ^(- a ^( at+H ^( i+ 5^ l (Hi) If - 1 



a, 



J k =k-i (** + 1 + + 1 + 3) (-«*k-, («*+*+ ?k-> (1 + |)i fc (I + Ik , 



„ (-««);„ (a„ +n+ + §),„_, (| + f ) 



i„it-, («» + » + ?k-i (i + fk(? + f ) 



2 "~' n'" 



(8) 



y(z) = LS a 



1 /a 



2 V2 



1 /a 



7 or- 



2 V2 



•1+7 



;M = ~P 2 ™ 2 (z,p);'? = -p 2 s« 4 fop) 



Wpi (-ab)/ (ab + |)/ S 




-«l).i("l + |).i(I).o(|)io ,-, 
-«l)/o(«l + 7)/o(?)ii(|)Ci 



2^7 (-a^),,(a A . + £ + |) ft (| + (1 + f ) ft . 



(-a*),,-! (<%+* + |Vi (I + |)%(1 + Ik 



(-«„),„ («„ + « + |k(f + !),„_. (1 + |k_i ■„ 1 



(9) 



([8]l is the first kind of independent solution of Lame function for the polynomial in 
Weierstrass's form as a — 2(2(Xj + j) or — 2(2a/ + j) — 1 where 7,0/ = 0, 1 , 2, • • • . 
And (O is the second kind of independent solution of Lame function in Weierstrass's 
form for the polynomial as a — 2(2a/ + j) + 1 or — 2(2a ; - + 7 + 1) where j, (Xj = 
0,1,2,-. 
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2.2 Infinite series 

The power series expansion of Lame function in algebraic form for the infinite series 
in Ref.0 is 



y(z) = ^Jn(z) =yo(x)+yi(x)+y 2 (x)+y3{x) + --- 





I I 




)i'o(4+2) ! 'o i'l=('o 1 ( 4 + 2 + 2 ) ! 'o ( 4 + 4 + 2 )<o(2 + 2 ) ! 'l (4 + 2 )< 

(2 fl -/7- c )(/ + |) 2 -^a(a + l) + ^(-f + |), (f +i+|) i0 



^a(a + l) + ^ 



(2a - Z? - c)(i k +i + |) 2 - £a(a + 1) + * 



(4 + 1+1+1)^ + 1+1+1) 

+ 2 + 2")i't( J + 2 + 4 + 2")'t(l + 2 + 2")'t-l (2 + 4" + 2")'<:-l 1 
4 T 2 T 2 "it-l v 4 ^ 2 ' 4 ' 2 /<*-l V 1 ' 2 ' 2 /<* V 2 ' 4 ' 2 J 
x £ (~f + | + |)'„(f + | + 4" + l)'^ 1 + | + + I + l n I 

i«=i n -l ( — 4 + I + 2") ! 'n-l (7 + f + 4 + 2")'/i-l + I + 2")'« (f + t + 2")'« J J 



Put © in (O. And put c = 1 as A=0 and co = (-p 2 ) 1,/4 as A = j into them. Then, 
we obtain two independent solutions of Lame equation in the Weierstrass's form. 

y(z) = LF (ju = -p 2 sn 2 (z,p);ri = -p 2 sn 4 (z,p)) 

°° r°- + iv 

_ ^ v 4/'oU t 4i | o ^i'n 
( =0 (l)'o(l)'o 

+ ^ £ -(i+p- 2 >f)+^(-f), (f + i), £ (-f+^(f + !),-, (j)^!)^ ,-, 

io=0 ( i 'o + 3)( / 0+4) U^oW'o ii=/ ( _ 7 + 2)io(f + 3)te(2)«i(3)»'i 

~ f - -(1+P~ 2 )'n + ^ f-SV C« + IV 

£ ~ (1 + p' 2) fa + fj! + gg (zi + |mi + k ly 1 + iv. (i + tVi ' 

>fc=i V ( l * + 1 + 1)(^ + 1 + 1) (-7 + 3)4-1 (7 + 3 + ( : + !)'*(! + §)>* , 
. f_« + 5v ra4.i-1.2v fi + »y (i 4. »y ' 

i n =t„-i \ 4 ' 2''n-l\4 ^ 4 T 2^/i-n 1 ' 2^'nU T 2- ,, « 



X 
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y(z) =LS(ii = -p 2 sn 2 (z,p);ri 



-p 2 STI 



°° (— _ 4- IV (_ 4- iV 



4 fep)) 



I 4 ' 4>to\ 4 T 2/'0 




L ( 4 + 4)11 ( 4 + (4)10(2 )'0 „li 



( 4 + 4)10 ( 4 + l)'o(4) £ 'l ( 2 )*1 



2 4 p 2 ( — f + 4")<o(? + 



2)* 



{4)10(^)10 



( 4 2 4 )jt ( 4 2 2 )'t ( 4 _ 2 )'*-! (1 + 2 )'*-! | 
( — f + 2 + 4)'*-l (f + 2 + 2") ! 't-l (4 + 2)'* (^ + 2)'*/ 

L (~T + 2 + 4)'"(t + 2 + 2)'"(4" + _)'"«—! (j + 2 )'"«-! 
f__4.54.IV f« + 2 + Iy f5.5Vfl4.5V 
i„=i„_! I 4^2^ 4^1-1 U T 2 T 2^"-l U ^ ^ 2< /i„ 



(12) 



(fT~) is the first kind of independent solution of Lame function in Weierstrass's form 
for the infinite series. And (IT_1 > is the second kind of independent solution of Lame 
function in Weierstrass's form for the infinite series. 
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3 Integral Formalism 



3.1 Polynomial in which makes B n term terminated 



The representation in the form of integral of Lame function in algebraic form for the 
polynomial in which makes B„ term terminated in Ref.[2| is 



y(z) = Ly»(z) 

n=0 

i =0 (1 + 2)10(4 + 2>'0 



-a ) 10 (a + i + A)i _ o 



n-l / r \ 



in 

n=l I i=0 



du„h u 1 



n-k u n -k 



(n-k-2+X) 



— I dv n - k — (l-V„_jt + i „V„-jfc(l-*„_fc)(l-«n-jfe)) ( " * +5+A) 

(v»_*-l) 1 



V n -jfc 1 - W„_Jt + i „V„_ J t(l — f„_^)(l — M„_fc) 



\{2a-b- c)V-^- k - 1+l) ( V„_, „«3y ) 2 

I V V n—k.p \ n K -. n vv / n—k,n 



(n-k-l+k) 



-a [ Ok-*-i + ~(n-fc-l+A) ) ( a„_ it _ 1 + -(n-fc--+A) ) +^ 



£° (-a ), (a + i + A) 1 - ( 



i'o=o (1 + 2)10(4 + 2 ) ( 'o 



(13) 



where 



f __1 V,- + ijv,-f,-M,- 

Vij = J (Vi- 1) l-Wljv ; (l-f,-)(l-«i) 
[ 77 only if i > j 



(14) 
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Put © in Cjj). 



y(z) = E^«W 

n=0 

1 



iff! (-ao), (ao + i + l),- ■ 

f ("fr 1 / , Un-k-l+X) f 1 , Un-k-2+X) 

+l i n /„ ^-*'»-t /„ <*»»-* 



= 1 I *=o 

] 

X 



70 



--; <fdv n ^ k -^— (1 -V„_ fc+ i ,„v„_ fc (l -*„_*) (1 -«„_*)) ( " * + * +A) 
27H J v„_* 



(v„-s- 1) 



Vn-jfe 1 — W„_t + i i „V„_ i t(l — f„_fc)(l - H„_t) 



(n-k-l+X) h 



2 4 P 



2 



i =0 U + 2V'oU+ 2i'0 J J 



Put c = 1 as A=0 and c = (-p 2 ) 1//4 as A = jj in ( TT3T >. 



1 /a 1 



jU = -p 2 sn 2 (z,p);ri = -p z 5 n 4 (z,p) 



2 „„4/ 



;'=o,i,2,- 



1 3 



2 F 1 [-a ,a + -;-;r 1 )+l t { ]J 



n-\ / ,1 



n=l I t=0 



dt n -k t n _ k 



o 



(n-k-2) 



*^ i jdv n -k : ^(l-^n-k+\,nVn-k{l-tn-k){l-^-k)) ( " 



(v n -t-l) 1 

Vn-k 1 — ^n-k+\.nVn-k{\ ~t n _ k )(l - U n _ k ) 



■II— k 



<{ - a +p-w:A ( r" 1) (^^) 2 **£ 



n-k-l) h 



2 4 p 2 



(16) 



( fTSI l is the integral formalism of the first kind of independent solution of Lame func- 
tion in Weierstrass's form for the polynomial as a = 2(20£,- + jf) or — 2(2a ; - — l 
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where j, a- s = 0, 1,2, • • • . 



1 ta 1 



y (z) = LS ttj ( a j = - [ - - - - j ) or - - ( - + 1 + j 



1 /a 



;jx = -p 2 sn 2 {z,p);r\ = -p l sn 4 (z,p) 



2 „,Ai 



2\2 

= V? <j 2 Fi (-ab, ao + I'^nj + £ I ft ( ^ fj 



7=0,1,2,- 



-A; 







Xr^ Idv„^ k -^— (1 - V„_ <:+ i i „v„_ fc (l -f„- J t)(l -«„_*)) ( " 



(v«-t-l) 1 

V«-jfc 1 — ^n-t+l,nV n -fc(l -?«-<:)(! ~ Un-k) 



n-k/i 24p2 



( fTvT > is the integral formalism of the second kind of independent solution of Lame 
function in Weierstrass's form for the polynomial as a = 2(2a ; - + l or — 2(2(Xj + 
7 + 1) where j, OCj = 0, 1,2, • • • . 



(17) 



3.2 Infinite series 

The representation in the form of integral of Lame function in algebraic form for the 
infinite series in Ref.j2] is 

y(z) = E?«( z ) 

n=0 

f 00 /'_ a j-Al /a i 1 i A ~\ 



i' =0 (1 + 2")'o(4 + l) 



'11 



»=1 I Ar=0 



^ <fdv„- k (l-V„_ fe+ i^V„_jfc(l-^_ ft )(l-M n _ J fe)) ( " 

27H J v n _ k 
(v n - k -l) 1 



Vn-k 1 - V'„_<; + l i „V n ^(l — f„-fc)(l — U n -k) , 



a n(n-^U-L q i I V ( f + l)'o(f + l + l)'Q ^'o 1 ,,n\ n on 
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Put © in (O. And put c = 1 as A=0 and co = (-p 2 ) 1,/4 as A = \ into them. Then, 
we obtain two independent solutions of Lame equation in the Weierstrass's form. 



y(x) = LF (fi = -p 2 sn 2 (z,p);ri = -p 2 sn A {z,p)) 

a a 13 \ ~ /-I Un-k-l) f 1 , l(n-k-2) 

■4'4 + 4 ; 4 ; T5(n(| dt ^Ck ^l^ k ut k 

— ldv n -i-^— (l-V„_fc+i ,„v n _jfe(l (l-u n - k )) ( " 
m J v„_ k 

(v„-*-l) 1 \i(f-+») 



X 

2ni 



v n-k 1 - ^n-k+l,nVn-k(l — ?«-*)( 1 — ««-*). 



yW = LS([i = -p 2 sn 2 (z,p);ri = -p 2 sn 4 (z,p)) 



1 



a 1 a 1 5 



n=l I k=0 



LI Tf I I' i ,(n-k-Z) f 1 , l{n-k-^ 

Sill /„ dt n-ktn-k J du n - k U n _ k 



X 



^ <fdv„- k -^— (l - tv > „_ i:+ i,„v„_i(l -f„_jt)(l -«„-*)) ( " * +?) 



x 



Vn-k 1 — V„_t + i „V„_i(l - f„_fc)(l — 



a la 15 



2*p 



2 



x ^(-fi'i+^ ;y ^ ^ (20) 



(fT9] l is the integral formalism of the first kind of independent solution of Lame func- 
tion in Weierstrass's form for the infinite series. And (120} is the integral formalism of 
the second kind of independent solution of Lame function in Weierstrass's form for 
the infinite series. 
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4 Asymptotic behavior of the function y(£, = sn 2 (z,p)) and the boundary 
condition for sn 2 (z,p) 

4. 1 Infinite series 



The condition of convergence in Lame function in the algebraic form for the infinite 
series and its asymptotic function in Ref. (2) are 

1 



limy(z) 

«>i 



1 



(2a— b—c) (x—a) 



(a—b)(a—c) ia—b)(a—c) 
(x — a) 2 (2a — b — c)(x — a) 



< 1 



(21) 



(22) 



(a — b)(a — c) (a — b)(a — c) 
Put (HI in (fJTJ and d22l , And its asymptotic function and the boundary condition of 
t, = sn 2 (z,p) for the infinite series of Lame function is 

1 



limv(,« feP))- 1+p2 , n 4 (Z)P) _ (1+p2)OT2(Z)P) 

where | (1 + p 2 )sn 2 (z,p) - p 2 sn 4 (z,p) \ < 1 

-2-- 2 



(23) 



-2 I 1+P 

More precisely, as < p < 



< sn 2 (z,p) < p 



-2 



(24) 



4.2 The case ofp «0 

Let assume that p is approximately close to 0. But p ^ 0. Then B n terms are negligi- 
ble. 

The condition of convergence in Lame function in the algebraic form for the case 
of 2a — b — c>-lor2a — b — c<§; — 1 and its asymptotic function in Ref. |2] are 

limy(z) = y -r- where a^b and c (25) 



1 



(2«— b— c)(x— a) 
(a—b)(a—c) 



(2a — b — c)(x — a) 



< 1 



(26) 



(a — b)(a — c) 

Put (HI in ( f25l l and d26l l, And its asymptotic function and the boundary condition of 
t, = sn 2 (z,p) for the polynomial in which makes B„ term terminated is 

1 



limy(sn 2 (z,p)) = 
The condition of convergence of sn 2 (z,p) is 



» 1' ' 1 - (1 +p 2 )sn 2 (z,p) 



< sn 2 (z,p) < 



1+p^ 



(27) 



(28) 
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5 Application 

Lame equation appears elsewhere in mathematical physics. For example, Recently, 
in "Droplet nucleation and domain wall motion in a bounded interval" |4), the au- 
thors investigate an extended model (a classical GinzburgLandau model) of noise- 
induced magnetization reversal. Lame equation arises in some specific boundary con- 
ditions. (see (8), (9) in Ref.QJ. In (9) its solution consists of the Jacobi eta, theta, 
and zeta functions according to Hermites solution of the Lame equation.) In "Group 
Theoretical Properties and Band Structure of the Lame Hamiltonian"|9], the authors 
represent a group theoretical analysis of the Lame equation, which is an example of 
a SGA band structure problem for su(2) and su(l, 1). (see (1), (10), (13), (14), (28), 
(29), (33), (38) in Ref.[9|) Applying three term recurrence formula|[T), we can obtain 
the power series expansion in closed forms and asymptotic behaviors of Lame func- 
tion analytically. And it might be possible to obtain specific eigenvalues for the Lame 
Hamiltonian. Again Lame equation is applicable to diverse areas such as theory of 
the stability analysis of static configurations in Josephson junctions [5 1, the computa- 
tion of the distanceredshift relation in inhomogeneous cosmologies 1 6 1 , magnetostatic 
problems in triaxial ellipsoids[7 | and etc. 



6 Conclusion 

From the above all, applying three term recurrence formula [Q], I show the power 
series expansion in closed forms of Lame function in the Weierstrass's form (infinite 
series and polynomials) and its integral forms. I show that a 2^1 function recurs in 
each of sub-integral forms of Lame function in the Weierstrass's form: the first sub- 
integral form contains zero term of A' n s, the second one contains one term of A„'s, the 
third one contains two terms of A„'s, etc. And I show asymptotic expansions of Lame 
function for infinite series and the special case as p ss 0. Since we obtain the closed 
integral forms of Lame function in the Weierstrass's form, Lame function is able to 
be transformed to other well-known special functions analytically; hypergeometric 
function, Mathieu function, Lame function, confluent forms of Heun function and 
etc. 

Various authors argue that the value of hp~ 2 can be chosen properly such that 
the Lame function is not an infinite series but a polynomial as a parameters of Lame 
functions is a positive integer. In the analysis of the three term recurrence formula |fl], 
since a is 2(2a ( - + i + X ) or — 2(2a, + i + X ) — 1 as i = 0, 1 , 2, ■ ■ ■ , Lame functions will 
be polynomial in which makes B n term terminated from d§) and ((9). X is the indicial 
roots which are or j, then all possible a is • • • , —3, —2,-1,0, 1,2,3, •• •. 

According to Erdelyi (1940[8 1), "there is no corresponding representation of sim- 
ple integral formalisms of the solutions in ordinary linear differential equations with 
four regular singularities; Heun equation, Lame equation and Mathieu equation. It ap- 
pears that the theory of integral equations connected with periodic solutions of Lame 
equation is not as complete as the corresponding theory of integral representations 
of, say, Legendre functions." The reason, why the analytic integral forms of Lame 
functions can not be obtained, is that the coefficients in a power series expansions 
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do not have two term recursion relations. We have a relation between three different 
coefficients. By using the three term recurrence formulaflT] , we are able to obtain an- 
alytic integral solution of any linear ordinary differential equation in which has three 
term recursion relations. 

7 Series "Special functions and three term recurrence formula (3TRF)" 

This paper is 7th out of 10. 

1 . "Approximative solution of the spin free Hamiltonian involving only scalar po- 
tential for the q — q system" lfl4l - In order to solve the spin-free Hamiltonian with 
light quark masses we are led to develop a totally new kind of special function the- 
ory in mathematics that generalize all existing theories of confluent hypergeometric 
types. We call it the Grand Confluent Hypergeometric Function. Our new solution 
produces previously unknown extra hidden quantum numbers relevant for descrip- 
tion of supersymmetry and for generating new mass formulas. 

2. "Generalization of the three-term recurrence formula and its applications" lfl5ll 

- Generalize three term recurrence formula in linear differential equation. Obtain the 
exact solution of the three term recurrence for polynomials and infinite series. 

3. "The analytic solution for the power series expansion of Heun function" lfl6ll 

- Apply three term recurrence formula to the power series exp ansion in closed forms 
of Heun function (infinite series and polynomials) including all higher terms of A„s. 

4. "Asymptotic behavior of Heun function and its integral formalism", ifTTl - Ap- 
ply three term recurrence formula, derive the integral formalism, and analyze the 
asymptotic behavior of Heun function (including all higher terms of A„s). 

5. "The power series expansion of Mathieu function and its integral formalism", 
IfTSI - Apply three term recurrence formula, analyze the power series expansion of 
Mathieu function and its integral forms. 

6. "Lame equation in the algebraic form" |[T9l - Applying three term recurrence 
formula, analyze the power series expansion of Lame function in the algebraic form 
and its integral forms. 

7. "Power series and integral forms of Lame equation in the Weierstrass's form" 
ll20l - Applying three term recurrence formula, derive the power series expansion of 
Lame function in the Weierstrass's form and its integral forms. 

8. "The generating functions of Lame equation in the Weierstrass's form" lETI - 
Derive the generating functions of Lame function in the Weierstrass's form (including 
all higher terms of A„'s). Apply integral forms of Lame functions in the Weierstrass's 
form. 

9. "Analytic solution for grand confluent hypergeometric function" 11221 - Apply 
three term recurrence formula, and formulate the exact analytic solution of grand 
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confluent hypergeometric function (including all higher terms of A„'s). Replacing /J. 
and eco by 1 and — q, transforms the grand confluent hypergeometric function into 
Biconfluent Heun function. 

10. "The integral formalism and the generating function of grand confluent hy- 
pergeometric function" ll23l - Apply three term recurrence formula, and construct 
an integral formalism and a generating function of grand confluent hypergeometric 
function (including all higher terms of A„'s). 
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